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. Abstract. The surprising lack of a high energy cutoff in the cosmic ray 

spectrum at the highest energies together with an apparently isotropic 
^ ■ distribution of arrival directions have strongly challenged most models 

proposed for the acceleration of ultra high energy cosmic rays. Young 
neutron star winds may be able to explain the mystery. We discuss this 
On I recent proposal after summarizing the observational challenge and plau- 

sible acceleration sites. Young neutrons star winds differ from alternative 
models in the predictions for composition, spectrum, and angular distri- 
' bution which will be tested in future experiments. 
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0\ . 1. Introduction 

ON 

The detection of cosmic rays with energies above 10^*^ eV has triggered consid- 
Qh' erable interest on the origin and nature of these particles. Many hundreds of 

Q . events with energies above 10^^ eV and over a dozen events above 10^'' eV have 

H I now been observed by a number of experiments such as AGASA (Hayashida et 

^ '. al. 1994, Takeda et al. 1998, Takeda et al. 1999), Fly's Eye (Bird et al 1993, 

1994, 1995), and Haverah Park (Lawrence, Reid, &: Watson 1991). Most unex- 
pected is the significant flux of events observed above 5 x 10^^ eV (Takeda et al. 
1998) with no sign of the Greisen-Zatsepin-Kuzmin (GZK) cutoff (Greisen 1966, 
Zatsepin & Kuzmin 1966). A cutoff should be present if the ultra-high energy 
\ particles are protons, nuclei, or photons from extragalactic sources. Cosmic ray 

protons of energies above 5 x 10^^ eV lose energy to photopion production off the 
cosmic microwave background and cannot originate further than about 50 Mpc 
away from Earth. Nuclei are photodisintegrated on shorter distances due to 
the infrared background (Puget, Stecker, & Bredekamp 1976) while the radio 
background constrains photons to originate from even closer systems. 

In addition to the presence of events past the GZK cutoff, there has been no 
clear counterparts identified in the arrival direction of the highest energy events. 
If these events are protons, they may point back to their sources within a few 
degrees, since at these high energies the Galactic and extra- galactic magnetic 
fields should not affect their orbits significantly. The gyroradius of a 10^*^ eV 
proton is 100 kpc in a /xGauss field which is typical for the Galactic disk, there- 
fore, protons propagate mainly in straight lines as they traverse the Galaxy. At 
present, no correlations between arrival directions and plausible optical coun- 
terparts such as sources in the Galactic plane, the Local Group, or the Local 
Supercluster have been found. The ultra high energy cosmic ray (UHECR) 
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data is consistent with an isotropic distributions of sources in contrast with the 
anisotropic distribution of Hght within 50 Mpc from Earth. 

2. The UHECR Puzzle 

In attempting to explain the origin of UHECRs, models confront a number of 
challenges. The extreme energy is the greatest challenge that models for astro- 
physical acceleration face, and to complete the puzzle, models have to match the 
spectral shape, the primary composition, and the arrival direction distribution 
of the observed events. 

If UHECRs are extragalactic, the observed highest energy event at 3 x 10^° 
eV (Bird 1994) argues for accelerators that reach as high as a ZcV (ZcV=10^^ 
eV) . The energetic requirements at the source increase with the distance traveled 
by the UHE primaries from source to Earth. Depending on the strength and 
structure of the magnetic field along the primary's path, the distance traveled 
may be significantly larger than the distance to the source. If 3 x lO^'^ eV is taken 
as a typical energy for protons travelling on straight lines, accelerators located 
further than 30 Mpc need to reach above 1 ZeV while those located further than 
60 Mpc require over 10 ZeV (Cronin 1992). As magnetic fields above ~ 10^^ G 
may thread intragalactic space (Riu, Kang &: Biermann 1998; Blasi, Buries &: 
Olinto 1999), protons travel in curved paths and sources need to be either more 
energetic or located closer to Earth. 

The extreme energy requirements have encouraged alternative explanations 
for UHECRs. For early universe physics, a ZeV is not particularly high in energy. 
For instance, relics from the Grand Unified scale at ~ 10^"^ eV may be causing 
these ultra high energy events (Hill 1983; Schramm & Hill 1983). The challenge 
for models that make use of early universe relics is generally the flux, the same 
challenge that observers face. At 10^^ cV, the observed flux of UHECRs is about 
~ 1 event/km^ /century which has strongly limited our ability to gather more 
than a dozen such events after decades of observations. Although challenging to 
observers and topological defect models, the flux is not particularly constraining 
in terms of general energetic requirements on astrophysical sources. In fact, this 
flux equals the flux of one gamma-ray burst that may have taken place in a 50 
Mpc radius volume around us (Waxman 1995; Vietri 1995). 

The energy spectrum of cosmic rays has a steep energy dependence N(E) oc 
E-'f, with 7 2.7 between ~ lO^eV and ^ lO^^eV and 7 ^ 3.1 for 10^^ 
^ E <^ 10^^ cV. Cosmic rays of energy up to ^ 10^^ eV are widely accepted 
to originate in shocks associated with galactic supernova remnants, but this 
mechanism has difficulties producing particles of higher energy. The events with 
energy above 10^^-^ eV, however, show a much flatter spectrum with 1 7 ^ 2. 
The drastic change in slope suggests the emergence of a new component of cosmic 
rays at ultra-high energies. This new component of cosmic rays is generally 
thought to be extragalactic, but they may also originate in an extended halo or 
in the Galaxy (Olinto, Blasi, & Epstein 1999) depending on their composition. 

The observed spectrum represents a convolution of the source primary spec- 
trum with the effect of energy losses during the propagation between the source 
and the Earth. If the new component is extragalactic, loss processes modify 
the spectral shape significantly. For primary protons the main loss processes 
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are pair production (Blumenthal 1970) and photopion production off the cosmic 
microwave background radiation (Greisen 1966; Zatzepin &; Kuzmin 1966). For 
straight hnc propagation, loss processes hmit sources of 10^^ eV to be within 
about 50 Mpc from us. For heavier nuclei, the infrared background induces losses 
(Stecker 1999) and the maximum distance for lO^'^ eV nuclei to originate from is 
~ 10 Mpc. In the case of photon primaries, the radio background represents the 
main source of losses and the distance is limited to ^ 10 Mpc. Depending on 
the poorly known strength and structure of the cxtragalactic magnetic fields, the 
GZK cutoff moves closer in distance for charged primaries. In addition, plausi- 
ble models have to accommodate the spectrum at the highest energies 10^° 
eV) without overproducing cosmic rays at ^ 10"^^ cV (Bcrczinsky, Grigorieva, &z, 
Dogiel 1990; Blasi & Olinto 1998; Sigl, Lemoine, & Biermann 1998). 

Charged particles of energies up to 10^° eV can be deflected significantly 
in cosmic magnetic fields. The Larmor radius of a particle with energy E and 
charge Ze in a magnetic field B, n = 0.1Mpc(£;/1020eV)/Z(5//xG). If the 
UHECR primaries are protons, only large scale intergalactic magnetic fields 
affect their propagation significantly, while for heavier nuclei the Galactic mag- 
netic field should also be taken into account. While the Galactic magnetic field 
is reasonably well studied, extragalactic fields are still very poorly understood. 

Together with a composition determination, the distribution of arrival di- 
rections can hold the key to the UHECR puzzle. Within a 50 Mpc radius volume 
around us, the most luminous structures are the Galactic plane, the Local Group 
and the general galaxy distribution with a relative overdensity aroTind the re- 
gion of the Local Supcrclustcr. If the UHECR source is dark matter, than the 
Galactic Halo is the relevant structure which is expected to be a spheroidal 
overdensity centered around the Galaxy. On larger scales the dark matter dis- 
tribution should correlate with the luminous matter distribution. For the few 
highest energy events, there is presently no strong evidence for correlations be- 
tween the events' arrival direction and any of these known structures (Stanev 
1999): the distribution is isotropic to first approximation. Since the number of 
observed events above lO^'^ eV is low, it is still early to tell. For slightly lower 
energies, some correlations have been detected. Recently, the AGASA group 
has announced that the distribution of 10^^ eV events shows a significant cor- 
relation with the Galactic Center and the nearby Galactic spiral arms (Yoshida 
et al 1999). If these correlations are confirmed, it would be strong evidence for 
a Galactic origin for cosmic rays around 10^^ eV. 

3. UHECR Accelerators 

As mentioned above, there are great difficulties with finding plausible accelera- 
tors for STich extremely energetic particles. Even the most powerful astrophysi- 
cal objects such as radio galaxies and active galactic nuclei can barely accelerate 
charged particles to energies as high as 10^*^ eV. 

Acceleration of UHECRs in astrophysical plasmas occurs when large-scale 
macroscopic motion, such as shocks and turbulent flows, is transferred to individ- 
ual particles. The maximum energy, £^max, is usually estimated by requiring that 
the gyro-radius of the particle be contained in the acceleration region. There- 
fore, -Emax is usually associated with the strength, B, and coherence length, L, of 
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Figure 1. B vs. L for E^aa; = P ZeB L 



the magnetic field embedded in the plasma, such that i?max ~ f3ZeB L, where 
usually (3 ^ v/c and Ze is the charge of the particle. As can be seen in Figure 
1, for Eiaax ^ 10^° eV and ZP ~ 1, the only known astrophysical sources with 
reasonable BL products are neutron stars {B ~ 10^^ Gauss and L ~ 10 km), 
active galactic nuclei (AGNs) {B ~ 10^ Gauss and L '-^ 10 AU), radio lobes 
of active galaxies {B ^ 10~^ Gauss and L 10 kpc), and clusters of galaxies 
{B ~ 10-6 Gauss and L ~ 100 kpc) (Hillas 1984). 

In general, when these sites are considered more carefully, one finds great 
difficulties due to either energy losses in the accelerating regions or the great 
distances of known sources from our Galaxy. In most of these objects one invokes 
shock acceleration as the primary acceleration mechanism. As discussed by 
Achterberg in this meeting, shock acceleration is not effective in reaching ZeV 
energies for most proposed accelerators with the possible exception of galaxy 
cluster shocks. The problem with clusters of galaxies as sources of UHECRs on 
Earth is that our location in the universe is not close enough to a cluster shock 
to avoid the GZK cutoff. Furthermore, UHECRs generated in typical clusters 
do not escape from them during the age of the universe (Blasi &; Olinto 1998). 

Moving left on Figure 1, radio lobes in FRII objects can reach the required 
energy if /3 ~ 1 (Rachen &; Biermann 1993). Again, location is a challenge since 
these are rare objects and far apart. A recent proposal, places the source of all 
UHECRs observed in a single object, M87, by invoking a Galactic wind that can 
make different observed arrival directions trace back to M87 (Ahn et al. 1999). 
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The existence of a Galactic wind with the required characteristics to allow for 
this possibility is yet to be determined. 

Active galactic nuclei (AGNs) are powerful engines as matter accretes onto 
very massive black holes. The problem with AGNs as UHECR sources is two- 
fold: one is the distance to more active objects and the other is common among 
highly energetic environments - losses due to intense radiation field downgrades 
particle energies well below the maximum achievable energy. These limitations 
have led to the proposal that remnant quasars, large black holes in centers of 
inactive galaxies, are UHECR accelerators (Boldt & Ghosh 1999). In this case 
one would have no obvious counterpart since any galaxy would be as likely 
to host such accelerators and losses are not as significant. The spectrum will 
be dominated by the local distribution of galaxies with more distance galaxies 
inducing a GZK cutoff (see, e.g., Medina-Tanco 1999). The detailed acceleration 
mechanism for this proposal is yet to be determined. 

Before discussing the last possibility in Figure 1, neutron stars, it is worth 
mentioning that the lack of a clear astrophysical solution for the UHECR puzzle 
has produced a number of models based on physics beyond the standard model 
such as monopoles, cosmic strings, cosmic necklaces, vortons, and superheavy 
long-lived decaying relic particles, to name a few. Due to the lack of space, we re- 
fer the interested reader to a few recent reviews (Berezinsky 1998, Bhattacharjee 
& Sigl 1999). 



3.1. Young Neutron Star Winds 

As shown in Figure 1, neutron stars may be effective in accelerating UHECRs 
(see, e.g., Berezinsky et al. 1990). Acceleration processes inside the neutron star 
light cylinder are bound to fail much like the AGN case: ambient magnetic and 
radiation fields induce significant losses (Venkatesan, Miller, &: Olinto 1997). 
However, the plasma that expands beyond the light cylinder is free from the 
main loss processes and may be accelerated to ultra high energies. 

One possible solution to the UHECR puzzle is our recent proposal that 
the early evolution of neutron stars may be responsible for the unexplained flux 
of cosmic rays beyond the GZK cutoff (Olinto, Epstein, & Blasi 1999; Blasi, 
Epstein, &; Olinto 1999). In this case, UHECRs originate in our Galaxy and 
are due to iron nuclei accelerated from the surface of strongly magnetic young 
neutron stars. 

Newly formed, rapidly rotating neutron stars may accelerate iron nuclei to 
ultra-high energies through relativistic MHD winds beyond their light cylinders 
(see, e.g., Michel 1991). The nature of the relativistic wind is not yet clear, 
but observations of the Crab Nebula indicate that most of the rotational energy 
emitted by the pulsar is converted into the flow kinetic energy of the particles 
in the wind. If most of the magnetic energy in the wind zone is converted into 
particle kinetic energy and the rest mass density of the wind is not dominated 
by electron-positron pairs, particles in the wind can reach a maximum energy of 
E,fnax — 8x10^'^ Z2QBisQ.^j^ eV, for iron nuclei {Z2Q = Z/26 = 1) and neutron star 
surface fields B = lO^^Bis Gauss and initial rotation frequency = SOOO^^k 
s~^. In the rest frame of the wind, the plasma is relatively cold while in the 
star's rest frame the plasma moves with Lorentz factors 7 ~ 

109 - iQio. 
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Figure 2. Solid lines for Ecr = 10^° eV and dashed lines for Ecr = 
3 X 10^° eV. The curves are plotted for two values of the envelope mass, 
Menv = 50 Msoiar and Menv = 5 M^oiar: as indicated. The horizontal 
line at spin period ~ 0.3 ms indicates the minimum period allowed for 
neutron stars. 



The typical energy of the accelerated cosmic rays can be estimated by con- 
sidering the magnetic energy per ion at the light cylinder Ecr — Bf^/STmcj 
where the Goldreich-Julian (1969) density is ncj = 1-7 x 10^^ iJiaOg^/Zcm"^. 
We find Ecr — 4 x 10^° ^26-613^^!,^. cV, similar to E^ax above. Therefore, neu- 
tron stars whose initial spin periods are shorter than ~ 4(55/10^^0) ms can 
accelerate iron nuclei to greater than lO'^^eV. 

About a year after the supernova explosion, the iron nuclei can escape 
through the remnant of the supernova that produced the neutron star without 
suffering significant spallation reactions. The supernova event ejects the enve- 
lope of the original star, making it possible for cosmic rays to escape. However, 
as the envelope expands, the young neutron star spins down and may become 
unable to emit particles of the necessary energy. A requirement for relativistic 
winds to supply UHECRs is that the column density of the envelope becomes 
transparent to UHECR iron before the spinning rate of the neutron star de- 
creases significantly. The allowed parameter space for this model is shown in 
Figure 2. 

The spectrum of accelerated UHECRs is determined by the evolution of 

the rotational frequency: As the star spins down, the energy of the cosmic ray 
particles ejected with the wind decreases. The predicted spectrum is very flat. 
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7 = 1, in agreement with the UHECR data. Furthermore, for the parameters 
within the allowed region, the acceleration and survival of UHECR iron nuclei 
is not significantly affected by the ambient photon radiation. 

Depending on the structure of the galactic magnetic field, the trajectories 
of the iron nuclei from galactic sources can be consistent with the observed 
arrival directions of the highest energy events (Zirakashvili et al. 1998). The 
gyroradius of these UHECRs in the Galactic field of strength Bg^i = 3/i Gauss 
is rs = 1.4 kpc£^2o/'^26 which is considerably less than the typical distance to 
a young neutron star 8 kpc). Furthermore, the cosmic ray component at 
10^® eV is nearly isotropic. If these cosmic rays are protons of Galactic origin, 
the isotropic distribution observed at these energies may be indicative of the 
diffusive effect of the Galactic and halo magnetic fields. Iron at lO^'^ eV probes 
similar trajectories to protons at a few times 10^^ eV. 



4. Conclusion 

Future experiments such as the Auger Project and the OWL-Airwatch satel- 
lite will be able to discriminate between different models (Cronin 1999; Watson 
1999). An excellent discriminator would be an unambiguous composition deter- 
mination of the primaries. In general. Galactic disk models invoke iron nuclei to 
be consistent with the isotropic distribution of events, and extragalactic astro- 
physical models tend to favor proton primaries, while photon primaries are more 
common for early universe relics. The observational tools in place for composi- 
tion discrimination are the muon content of shower in the ground arrays (more 
muons for nuclei vs. nucleons) and the depth of shower maximum in fluorescence 
detectors (the heavier the primary the deeper in the atmosphere their shower 
maximum). In addition, the correlation of arrival directions for events with en- 
ergies above 10^*^ eV with some known structure would be key in differentiating 
between different models. For example, we should see a correlation with the 
Galactic center and disk for the case of young neutron star winds (see, Stanev 
1999; Yoshida et al. 1999), and the large scale galaxy distribution for the case 
of quasar remnants. Both aspects will be testable with future experiments. 
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